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Low-temperature thermodynamic properties of the spin-i Heisenberg ladder system with non- 
magnetic impurities are discussed using an effective Hamiltonian for the impurity-induced spins 
in the background of a spin liquid with a gap. It is shown that the uniform susceptibility 
shows Curie-like behavior in high-, intermediate-, and low-temperature regimes, with three dif- 
ferent Curie constants. The ratio of specific heat to temperature, C/T, is divergent in the 
low-temperature limit, reflecting a singular distribution of the effective couplings among the 
impurity spins. 
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Quasi-one-dimensional quantum spin liquids with a 
spin excitation gap have attracted much interest for more 
than a decade. One recent example of interest is the spin- 
h ladder system because of its relationship to cuprate su- 
perconductors, jjj These spin ladders can be modeled by 
a nearest-neighbor Heisenberg Hamiltonian of the form 



|U=1,2 



i,2 



(1) 



for the case of a two-leg ladder, where Sj^ denotes the 
spin operator on rung i and leg [i (see Fig. 1(a)). Consid- 
erable effort has been invested both experimentally and 
theoretically to study the properties of these systems. 
We may say that by now, regular ladders, in particular, 
those with two legs, are well understood both qualita- 
tively and quantitatively. In general ladder systems with 
an even number of legs have a resonating valence bond 
(RVB) ground state with a spin gap, while an odd num- 
ber of legs leads to gapless excitations in a ground state 
with quasi-long-range order. The spin gap is largest for 
the system with two legs, where numerical data show 
A 0.5 J. H] Therefore, in real two-leg ladder systems, 
A can be very large, as the examples (VO)2P207 with 
A - 40K and SrCu 2 3 with - 400K show. § 

An interesting new aspect arises, if we take disorder 
into account. The type of disorder we are interested in 
here, is the depletion of spins. This may occur when 
ions delivering the spin degrees of freedom are replaced 
at random by other ions without spins, by accident or by 
intentional doping [Fig. 1(b)]. For example, in SrCu2 03 
some Cu-ions can be substituted by nonmagnetic Zn- 
ions. [| It has been observed experimentally that in 
systems with a dimer ground state, such a depletion in- 
troduces a strong component of staggered moments. Q 
This effect has recently been investigated theoretically. 
J^, ^, In the vicinity of a nonmagnetic impurity the 
spins develop a strong staggered correlation, generating 

* On leave from Theoretische Physik, ETH-H6nggerberg, 8093 
Zurich, Switzerland. 



an effective spin doublet. These effective spins dominate 
the low-temperature physics, since, due to the presence 
of the spin gap, all other spin degrees of freedom dis- 
appear at small enough energy scales. For a low, but 
finite concentration of impurities, the correlation among 
these effective spins will play an important role in the 
low-energy properties. It is the aim of this letter to 
show some consequences of this correlation which may 
be tested experimentally. 
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Fig. 1. (a) Regular spin ladder system, (b) depleted spin ladder 
system. 



Obviously each empty site (nonmagnetic impurity) in 
the RVB ground state of a two-leg spin ladder represents 
a spin £> = @.|7| In the following we will assume that 
the impurity concentration is sufficiently low such that 
the idea of a spin-i degree of freedom associated sepa- 
rately with each impurity is valid. The coupling between 
these impurity spins is very weak and short-ranged, be- 
cause it is mediated by the spins which participate in 
the formation of the RVB spin liquid state. Before ad- 
dressing the issue of coupling strengths we discuss a more 
qualitative aspect. It is intuitively clear that the coupling 
between the effective spins can be either ferromagnetic or 
antiferromagnetic, depending on their relative locations. 
We can show this on a rigorous basis by considering a 
ladder system with two impurities at sites (ii,//i) and 
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(«2,A*2)- The nature of the effective coupling can be de- 
termined by finding the total spin quantum number of 
the ground state. This can be straightforwardly done 
using Marshall's theorem || which states that the spin 
quantum number of the ground state of Ti. in eq. (Q) is 
identical to that of the Hamiltonian, 



H 



EE* 



(2) 



where A and B denote the sets of sites belonging to the 

two different sublattices (note that the ladder system is 

a bipartite lattice). The eigenenergies E of Ti. depend 

on three parameters: S to t (Stot = #A (Sa = 

J2ieA S i) and S b (S B = J2ieB S *)- Tne y are simply 
given by 

E(S t ot,SA, Sb) — g [Stot {Stot + 1) — Sa(Sa + 1) 

-S B (S B + 1)]. (3) 

The ground state corresponds to the case when Sa and 
Sb have their maximal and Stot its minimal value (= 
\Sa — Sb\)- This immediately leads to the well-known 
result that if the number of sites on the two sublattices is 
identical, then Sa — Sb and the ground state is a total 
spin singlet. 

We can now apply this theorem to our problem. If the 
two impurities are on different sublattices, they remove 
one site from each sublattice. Thus the ground state 
again has Sa — Sb and is a singlet, implying that the 
coupling between the two impurity spins is antiferromag- 
netic. On the other hand, if the two impurity sites lie 
on the same sublattice, then the spin quantum number 
of this sublattice is smaller by 1 than that of the other 
and Stot = 1, a triplet ground state. This corresponds to 
a ferromagnetic coupling. Consequently, the low-energy 
properties of the randomly depleted spin ladder may be 
described effectively by a model of spins (S = |) with 
couplings of random strength and random sign. The low- 
temperature properties of this type of system have re- 
cently been investigated by means of a real-space renor- 
malization group scheme |9|, [10| based on earlier studies 
by Ma and co-workers. |ll| , |l2| In the following we discuss 
how these results can be applied to this system. 

Let us first consider the uniform susceptibility of this 
system. At high temperatures (k^T > J) the spins are 
essentially free and the susceptibility per site obeys the 
Curie law, 



X(T) = 



Mb 



(4) 



4/cbT 

where /ib is the Bohr magneton and &b is the Boltzmann 
constant. It is clear that for temperatures below the 
spin gap, only the effective impurity spins contribute to 
the susceptibility, because the other spins have frozen 
out. This shows up as a sharp drop in the susceptibility 
for T < A. Since the impurity spins are only weakly 
coupled, they will yield a Curie-like susceptibility in an 
intermediate temperature regime: 



X(T) 



4fc B T' 



(5) 



where z is the impurity density (z < 1). At lower tem- 
peratures, these impurity spins start to correlate due to 
their interaction. First the two spins with the strongest 
mutual coupling freeze into either a singlet or a triplet 
state, depending on the sign of the interaction. In the 
case of the triplet configuration, they act as a single new 
spin (S = 1) coupled via weaker effective interactions 
to its two adjacent impurity spins. For the singlet case, 
however, these spins no longer contribute as degrees of 
freedom, but still mediate an effective interaction be- 
tween the two adjacent impurity spins through virtual 
excitations. In this way we find the equivalent effective 
spin system where the strongest bond of the original sys- 
tem was integrated out. With decreasing temperature, 
the next strongest pair forms a single spin with new effec- 
tive coupling to its neighbors and so on. By iteration of 
this process an effective spin system gradually evolves, 
in which spin degrees of freedom consist of large clus- 
ters of randomly correlated impurity spins. The forma- 
tion of these effective spins is determined by the energy 
scale &bT; all degrees of freedom with correlation en- 
ergies larger than this energy scale are frozen into such 
effective spins. Consequently these effective spins behave 
essentially independently, since the correlation energies 
among them are smaller than the thermal energy. Using 
this picture we can calculate the uniform susceptibility 
in the following way. 

Let us assume that the average size of the cluster of 
correlated spins contains n rungs. Note that n is a mono- 
tonic function of the temperature. Now we can estimate 
the average size of the effective spin associated with the 
cluster by again applying Marshall's theorem. Consider 
a Heisenberg ladder containing n rungs where 2zn spins 
are randomly depleted on average. The effective spin size 
of the cluster corresponds to the ground-state spin quan- 
tum number Stot of this finite-length Heisenberg ladder, 
which is, on average, 

(S M ) = (\S A -S B \) = ±(\N A -N B \), (6) 

where Na (Nb) is the number of depleted sites on the 
A (B) sublattice and the average is taken over the im- 
purity configuration. This can be easily estimated from 
a random- walk picture: 



{Stot) - vf 



(7) 



for nz 1. Thus the uniform susceptibility per site 
becomes 



X(T) = 



A (SL(T)) 



(8) 
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This result means that the susceptibility also follows a 
Curie law at sufficiently low temperatures. However, the 
Curie constant for T — > is different from that in eq. (||). 
Consequently, we expect to see crossovers in the Curie 
constant, as schematically shown in Fig. 2. Our predic- 
tion is that the Curie constant drops by a factor of 1/3 
from the intermediate- to the low-temperature regime. 

The characteristic temperature T* where this 
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Fig. 2. Schematic picture of the three regimes of the Curie con- 
stant, yT. The first crossover temperature is the spin gap A and 
the second, T* , is related to the average coupling strength. 



crossover occurs is set by the average coupling strength 
among the impurity spins (Fig. 2). Here we attempt to 
give a rough estimate of this energy scale based on a 
simplified model of the spin gap state. Suppose only two 
spins, at sites (0, 0) and (j, fi), are depleted in the Heisen- 
berg ladder, and the resulting impurity spins interact 
with each other by exchanging a spin triplet excitation 
in the RVB spin liquid. For calculating the long-distance 
behavior of this interaction, it is sufficient to consider an 
approximate dispersion of these excitations. We adopt 
the dispersion given by Troyer et al., which they used 
to fit numerical data of the regular two-leg spin ladder 
system. [fl3f Their momentum dependence for the spin 
triplet excitation energy has the form 



E k = ^A 2 +V 2 (k 



(9) 



where k is the momentum along the leg. Note that this 
excitation corresponds to an odd-parity state along the 
rung. The effective interaction can be obtained, in a 
similar way to the RKKY treatment, by second-order 
perturbation: 

n J 2 p ik i 
JU,») « — (-l) J+Al / dfc - 

" ^ J- n ^(A/v) 2 +k 2 

(-iy+»K (j/i ), (io) 



v 

2aJ 2 



where a is an unknown constant of order unity com- 
ing from the matrix elements. The parameters used by 
Troyer et al. are A « 0.5 J and v = 1.8 J. |1| Thus we 
obtain 



J(j^)^l.laJ(~iy+^K U/£ ). 



(11) 



It is important that the length scale appearing in the 
modified Bessel function Kq, Iq — v/A « 3.6, is rather 
short, and equivalent to the spin-spin correlation length 
in the regular ladder. |l3| Therefore it seems reason- 
able to assume that the effective spin system for low- 
temperature physics (T < A) is described well by con- 
sidering only nearest-neighbor Heisenberg coupling. 

For a given z the distribution of distance between 
neighboring impurities is V{t) — 2zexp(— 2z£), which 
gives an average distance {£) of l/2z. This leads to an 



average coupling strength of 



(|J|) m l.laJ,/^zT Q e' 1/2zl0 . 



(12) 



- 5-10" 2 aA, 
400K, as in 



Thus, for 5% of depleted sites we get (| J|) 
which could be on the order of 10K for A 
SrCu 2 3 . 

So far we have ignored the fact that ladder sys- 
tems can be easily disconnected if we remove sites ran- 
domly. Therefore the above analysis is valid only when 
a short-range correlation among the impurity positions 
is present, which prohibits configurations leading to de- 
coupled segments of the ladder. If we ignore this type of 
correlation, then the ladder decays into decoupled seg- 
ments of finite length, which, on average, contain 1/z 2 
rungs and 2/z impurities. This average length l/z 2 
serves as a low-energy cutoff. This means that the pro- 
cess of developing spin correlation with decreasing en- 
ergy scale is interrupted once all spins in a segment are 
frozen into its ground-state configuration. Hence in the 
low-temperature limit, each segment behaves as an in- 
dependent effective spin degree of freedom, yielding a 
Curie-like temperature dependence of the susceptibility. 
The Curie constant is naturally expected to be some- 
where between the values given in eqs. §) and (§). We 
can estimate this Curie constant by calculating the mean 
value of the ground-state spin quantum number S' of the 
segments. For a given number n of impurity spins in a 
segment, we obtain 

1/2 



S'(n) 




1/2 



(13) 

where Sf denotes impurity spins ±1/2. Noting that the 



distribution of n is given by Vin) 
obtain the Curie constant 

^ B (S'(n)[S'(n) 
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(14) 



which gives only a small deviation from eq. (|J) for z«l. 

Next we discuss the low-temperature behavior of the 
specific heat. For this purpose we need to know the dis- 
tribution of the effective couplings J. The fact that the 
coupling strength decays exponentially with the distance 
leads to a distribution which is rather singular for J — > 0. 
In this limit, it is approximated by 



d£V(£)8(\J\- JoVV^e-^") 

2z£ -l 




(15) 



where we have neglected small logarithmic corrections. 
The smaller the concentration z, the more of a singular 
distribution we obtain. It was shown by Westerberg and 
co-workers || that in a real-space renormalization group 
treatment, the distribution of couplings approaches a 
universal singular form E3| for J — > 0. This leads to 
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a universal low-temperature behavior of various quanti- 
ties such as the specific heat; at very low temperatures, 
the specific heat C would follow a power law, and the 
ratio C(T)/T would be more singular than that of a fer- 
romagnet (C fcrro /T oc T" 1 / 2 for T — ► 0). This can be 
easily seen from the following argument. 

As we have already discussed, at low temperatures 
many impurity spins correlate in clusters which then be- 
have more or less as single spin degrees of freedom. The 
number of correlated impurity spins, n(^> 1), at tem- 
perature T can be estimated by noting that the finite- 
size gap, An, to the first excited state in this n spin 
system should be on the order of k^T. It is natural 
to assume that A„ oc n -1 /", where the exponent a 

From eq. (Q) we 



get S'(n) oc A„ Q ^ 2 oc can then estimate 



should be determined from X>(|J|). 

oc 

the entropy a per site by assuming that each effective 
spin is essentially independent: a = zln[2S"(n) + l]/n. 
From these relations, we get cr{T) oc zT a \n(l/T), yield- 
ing C(T) oc zT Q ln(l/T) to leading order in T. If the 
couplings J between the impurity spins were all uniform 
and ferromagnetic, C(T) oc T 1 / 2 . With random distri- 
bution of coupling, we expect that the low-lying excita- 
tions are, in general, softer, because they are dominated 
by spin fluctuations around the weakest couplings. This 
leads to the conclusion that a < 1/2, in agreement with 
ref. I Hence, C{T)/T > T" 1 / 2 for T -> 0. 

From these arguments we can expect the following 
temperature dependence of the specific heat. At high 
temperatures the specific heat behaves essentially like 
that of the regular Heisenberg ladder. It has a peak at 
k^T ~ J and drops towards lower temperature due to 
the presence of the spin gap. Note that the staggered 
spin environment of each impurity introduces a discrete 
spectrum of localized excitations in the uniform spin gap, 
which can smear the gap behavior. As the tempera- 
ture is decreased further, it shows a small broad peak 
or shoulderlike structure around k^,T w (| J|) due to the 
interaction between the impurity spins, and it becomes 
the anomalous power law behavior discussed above for 
T — > 0. However, in the case of uncorrelated impurity 
positions, this power law behavior will be cut off by the 
finite-size gap determined by the average length 1/ z 2 and 
the coupling distribution 2?(J). Obviously the ground 
state of the whole system will have a large degeneracy, 
unless there is residual interaction between "decoupled" 
segments. 

In summary, we have found that at low temperatures 
the randomly depleted spin ladder system behaves in 
many respects like a random spin system with ferro- 
and antifcrromagnetic couplings. Nevertheless, we em- 
phasize that the impurities induce a staggered correla- 
tion which is, although inhomogeneous, coherent on a 
long length scale at low temperature. This would be the 
dominant spin correlation feature in neutron scattering 
measurements. An important result of our study is that 
in this system correlation effects lead to three distinct 
temperature regimes. Of particular interest is the low- 
temperature regime, where Curie behavior of the uniform 
susceptibility appears with a nontrivial Curie constant. 



Furthermore, the specific heat shows an anomalous tem- 
perature dependence for T — ► 0, as long as decoupling 
into segments is not significant. These properties are 
not restricted to two-leg ladders, but would also appear 
in any even-leg ladder and other systems with an RVB 
ground state. We believe that for a certain range of 
impurity doping in ladder systems, these effects would 
be in an experimentally accessible temperature range, in 
particular, for two-leg ladders with a large spin gap. In 
actual Heisenberg ladder systems such as SrCu203, there 
is certainly coupling between different ladders. Although 
very weak, this coupling can also introduce a low-energy 
cutoff for the temperature dependence of \ and C and 
can lead to long-range order of the staggered spin cor- 
relation, as was actually observed experimentally. 
This change to 3-dimensional behavior would spoil the 
observation of our low-temperature regime. 

Acknowledgements 

The authors thank N. Nagaosa, T. M. Rice and B. 
Frischmuth for helpful discussions, and H. Fukuyama for 
providing a preprint of his paper prior to publication. 
They are also grateful for financial support by Monbusho 
which made this collaboration possible. One of the au- 
thors (M.S.) acknowledges a fellowship (PROFIL) by the 
Swiss Nationalfonds. 



[l] 



E. Dagotto and T. M. Rice: Science 
Azuma, Y. Fujishiro and M. 



agi, M 



See for a recent review 
271 (1996) 618. 
[2] M. Nohara, H. Takagi 

Takano: preprint. 
[3] M. Azuma, Y. Fujishiro, M. Takano, T. Ishida, K. Okuda, M 
Nohara and H. Takagi: preprint. 

L. P. Regnault, J. P. Renard, G. Dhalenne and A 
Revcolevschi: Europhys. Lett. 32 (1995) 579. 
H. Fukuyama, T. Tanimoto and M. Saito: J. Phys. Soc. Jpn 
65 (1996) 1183. 

H. Fukuyama, N. Nagaosa, M. Saito and T. Tanimoto: per- 



[1] 
[5] 
[6] 



print. 
[71 G. B. 



Mart ins. E. Dagotto, and J. A. Riera: 



mat/9605069 



:ond- 



i] See, for example, A. Auerbach: Interacting Electrons and 
Quantum Magnetism (Springer- Verlag, New York, Berlin, 
Heidelberg, 1994). 
[9] E. Westerberg, A. Furusaki, M. Sigrist and P. A. Lee: Phys. 
Rev. Lett. 75 (1995) 4302. 
[10] See also A. Furusaki, M. Sigrist, E. Westerberg, P. A. Lee, K. 

B. Tanaka and N. Nagaosa: Phys. Rev. B 52 (1995) 15930. 
[11] S.-K. Ma, C. Dasgupta and C.-K. Hu: Phys. Rev. Lett. 43 
(1979) 1434. 

[12] C. Dasgupta and S.-K. Ma: Phys. Rev. B 22 (1980) 1305. 
[13] M. Troyer, H. Tsunetsugu and D. Wiirtz: Phys. Rev. B 50 
(1994) 13515. 

[14] For high degree of singularity in the limit | J| — > 0, however, 
we generally cannot expect that T>(\ J\) approaches a universal 
scaling regime. According to eq. (15) small values of z may 
create such situations. 

[15] A recent example may be CaV^Og, which shows a plaquette 
RVB ground state; see K. Ueda, H. Kontani, M. Sigrist and 
P. A. Lee: Phys. Rev. Lett. 76 (1996) 1932. 



